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Abstract 

A qualitative analysis is presented for a class of homogeneous cosmologies de- 
rived from the string effective action when a cosmological constant is present 
in the matter sector of the theory. Such a term has significant effects on the 
qualitative dynamics. For example, models exist which undergo a series of os- 
cillations between expanding and contracting phases due to the existence of 
a heteroclinic cycle in the phase space. Particular analytical solutions corre- 
sponding to the equilibrium points are also found. 
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I Introduction 



Very early universe cosmology provides one of the few environments where the pre- 
dictions of fundamental theories of physics, and in particular string theories, can be 
investigated. String theory is the most promising candidate for a unified theory of the 
fundamental interactions. It introduces significant modifications to the standard, hot 
big bang model based on conventional Einstein gravity and a study of string-inspired 
cosmologies is therefore important. 

String theories predict the existence of a graviton, f/^,^, a scalar 'dilaton' field, 
and an antisymmetric two-form potential, Bf^,^, with a field strength Hfj^i^x = d[^B^x] 

0]. In four dimensions, the three-form field strength is dual to a one-form, V^a, 
such that H'^'^^ = e'^e^'^^'^V ^a, where ^^'^^'^ is the covariantly constant four-form [Q. 
The one-form may be interpreted as the gradient of a scalar 'axion' field. The string 
field equations can then be derived from the effective action 



5* = y d^x^J—ge 



/2+(V$)2-ie2* (Va) 
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M, 



where Sm represents the action for perfect fluid matter sources, R is the Ricci cur- 
vature of the spacetime and g = dett^^,^. The dilaton-graviton sector of action (||) 
may be interpreted as a Brans-Dicke theory, where the coupling parameter between 
the two fields takes the specific value u; = — 1 [Q. The value of the dilaton field 
determines the effective value of Newton's 'constant', Ges oc e*. 

The general solutions to the field equations of action ([^) are known analytically 
when Syi = for both the spatially flat and isotropic Friedmann-Robertson- Walker 
(FRW) universes and the anisotropic Bianchi type I models ^. The purpose of 
the present paper is to qualitatively investigate the consequences of introducing a 
cosmological constant. Am, into the matter sector of Eq. (|1|): 

S = J d'x^ |e-* [i? + (V*)' - ie^* (Va)^] - Am} . (2) 

This term may be interpreted as a perfect fluid matter stress with an equation of state 
p = —p. It could be generated by a slowly moving scalar field, with a kinetic energy 
contribution dominated by a self-interaction potential, p ~ ~V ~ —p. Analytical 
FRW solutions have not been found for this model when the axion field is trivial and 
Am > 0. Moreover, the combined effects of the cosmological constant and axion 
field have not been considered previously. 

We determine the general structure of the phase space of solutions for spatially 
flat FRW and axisymmetric Bianchi type I cosmologies derived from action @ for 
arbitrary Am. This complements the work of Refs. 0, |rT|, [12], |T3|, where the quali- 
tative effects of introducing a cosmological constant. Am oc e~*, into the gravitational 
sector of Eq. ([^) were determined. 

The paper is organized as follows. In Section 2, the cosmological field equations 
and solutions for a zero cosmological constant are presented. The qualitative be- 
haviour of the models with positive and negative Am is determined in Sections 3 and 
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4, respectively. The phase portraits are interpreted in Section 5 and we conclude with 
a discussion in Section 6. 



II Cosmological Field Equations 

The metric for the Bianchi type I model may be written in the form 

ds^ = -df + habdx''dx\ a,b = l, 2, 3, (3) 

where hab{t) is a function of cosmic time t only and represents the metric on the 
surfaces of homogeneity. The axisymmetric model may be parametrized by hab = 
g2a(i) |^g2/3(t)^ where e'^" denotes the effective spatial volume of the universe. The 
traceless, diagonal matrix Pat = diag[/?, /?, — 2/?] determines the shear of the models 
and we refer to (3 as the shear parameter ||I4|. The spatially flat, isotropic FRW 
model is recovered in the limit where (3 = and, in this case, represents the scale 
factor of the universe. 

Substituting the metric (H) into the action (H) and integrating over the spatial 
variables implies that 

S = Jdt e^" |e~* 6d<i> - 6a^ + - + ^e^^a^ - Am} , (4) 

where the co-moving volume has been normalized to unity without loss of generality 
and a dot denotes differentiation with respect to t. The field equations derived from 
Eq. are given by 

a = d0 + 0' - - 6/3^ - ^AMe^+=^° (5) 

^ = + 6/^2 + lAMe'^+=^" (6) 

a = —{ip + 6a)& (7) 
(3 = (8) 

= $ - 3a (9) 



where 



defines the 'shifted' dilaton field and the generalized Friedmann constraint takes the 
form 

3«2 - 02 ^ 6/52 + i(T2e2^+6" + AMe^+='" = 0. (10) 
Eqs. (|5|)- (p!oD may be simplified by introducing the new time coordinate 

-(v+3a)/2^ 

de- dt ^ ' 
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and employing the generalized Friedmann constraint equation (|TOD to eliminate the 
axion field. The remaining field equations are then given by 



a 

a" 



1 



(3' - 3a') 



(12) 
(13) 
(14) 



where a prime denotes differentiation with respect to 6. 



The general solution to Eqs. (^-(|T3) is known when the cosmological constant 
vanishes 0. It is given by 



1/2 
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\s/s^\ + |s/s* 



(15) 



where s = /* cit'e""^*'^ is conformal time, {a^,, s^,,^^,a^,, /3^,} are arbitrary constants 
and {r, q} satisfy the constraint equation r = (3 — 12g^)^/^. 

The solutions to Eqs. (^^(0) for a trivial axion field and zero cosmological 
constant have a power-law form: 



e \t\ 

<S>, I ,|±3/i.-l 



,n±y/{l-3hi)/6 



(16) 



where h^: is a constant such that < ly/S. The solution (|I5|) asymptotes to these 
power-law models at early and late times and the axion field is therefore dynamically 
negligible in these limits. When an axion field it present, as in Eq. ([I5|) , the universe 
undergoes a smooth transition between the two power-law solutions (|T^) and exhibits 
a bounce when s ~ s^,. In the isotropic limit, hi = 1/3, and the time-reversal of the 
e° oc solution is infiationary. It corresponds to the pre-big bang cosmology, 

where the infiationary expansion is driven by the kinetic energy of the dilaton field 

In the next section we determine the phase portraits for the generalized model 
with a non-trivial axion field and Am > 0. The effect of the cosmological constant 
on the solutions (pTSf) can then be established. 
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Ill Positive Cosmological Constant 



When Am > 0, we can rewrite Eqs. ([T^)-(p^ using new variables defined by 

h = a\ ^ = ^\ N = I3'. (17) 

Eq. (p!OD then imphes that 

^/-^ > + 6A^^ + Am > (18) 
and consequently we may normalize with if). We therefore define 

X = —— 19 

y - ^ (20) 

- - -- (22) 

and assume that ^/^ > 0. (The case '0 < is related to a time-reversal of the system 
and the qualitative behaviour is similar). The three-dimensional system (|T^)-(|T^) is 
therefore given by 

^ = (a; + v/3)[i_a;2_^_^] + l^[a;_y3] (23) 



^ = 2i/|[l-x^-y-;.] + -2;| (24) 
^ = 2z[[l-x'-y-z]-\{l-z-V3x)]. (25) 

It follows from the definitions (p!9[)-(pTD that the phase space is bounded with 

< {x'^,y,z} < 1 subject to the constraint 1 — x"^— y — z>0. The invariant set 

1 — x'^—y — z = corresponds to a zero axion field. The dynamics of the system 



(23)-(25) is determined primarily by the dynamics in the invariant sets y = and 
z = 0. These correspond to a zero shear parameter and a zero cosmological constant, 
respectively. The dynamics is also determined by the fact that the right-hand side of 
Eq. (|2^) is positive-definite so that y is a. monotonically increasing function. This 
guarantees that there are no closed or recurrent orbits in the three-dimensional phase 
space. 
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Ilia Isotropic Model for Am > 

The isotropic FRW cosmology corresponds to the invariant set y = 0, where the shear 
parameter is trivial The system (P5D-(P5D reduces to the following plane system in 
this case: 

^ = {x + V3)[l-x'' -z] + ^z[x-V3] (26) 
^ = 2z[[l-x'-z]-^-{l-z-V3x)Y (27) 

The equilibrium points and their associated eigenvalues are given by 

S,: x = -l,z = 0; Ai = 2(^3-1), A2 = -(1 + Vs) (28) 
x = l,z = 0; Ai = -2(V3 + 1), A2 = (V3-1) (29) 

^- ^ = -^'^=i; Ai,2 = ^±^v^. (30) 

The points 5*1 and S2 are saddles and F is a repelling focus. The phase portrait is 
given in Fig. 1. 

z 




Figure 1: Phase portrait of the system l^^)-i^2li), corresponding to the isotropic FRW 
model with Am > 0. Equilibrium points are denoted by dots and the labels in all 
figures correspond to those equilibrium points (and hence the exact solutions they 
represent) discussed in the text. We shall adopt the convention throughout that large 
black dots represent sources (i.e., repellors), large grey- filled dots represent sinks (i.e., 
attractors), and small black dots represent saddles. Note that in this phase space orbits 
are future asymptotic to a heteroclinic cycle. 

In the invariant set 1 — — 2; = 0, corresponding to the case of a zero axion field, 
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Eqs. ( p6D and ( j27|) reduce to the single ordinary differential equation: 

§ = i(l-.^)(.-V3). 

which can be integrated to yield an exact solution in terms of B-time. 



(31) 



Illb Anisotropic Model for Am > 



In the full system (P3D-(P5D, corresponding to the anisotropic model with a non-trivial 
shear parameter, there exists the isolated equilibrium point (and their associated 
eigenvalues) 



1 



X 



W 



3v^' 
(Ai, Ai, A3) 
?/ = 1 - x^ 
(Ai, Ai, A3) 
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(32) 
(33) 



Hence, F is a global source. The set W lies in the invariant set ^ = on the boundary 
y = 1 — x . Points on W with x E (-1/v^, l/VS) are local sinks, while the remaining 
points are saddles in the full three-dimensional phase space. (In the invariant set 2; = 
equilibrium points with x E [—1, — 1/a/3) are repelling and those with x E {—l/y/S, 1] 
are attracting). The phase portrait for this system is given in Fig. 2 and table |l| lists 
each equilibrium set and its stability.. 



Equilibrium Point 


Stability 


F: Eq. (H) 


Repellor (Source) 


W: Eq. (ID 


Attractor (Sink) for ^ < x < 
Saddle otherwise 



Table 1: Equilibrium sets for anisotropic model with Am > 0, and their stability (the 
equations where each set is defined is also listed). 



We note that there exists an exact, anisotropic solution of Eqs. 

1 



X 



3^3 



constant 



and 
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y = — ^ I ^ 



16 
27 



-(psD, where 
(34) 

(35) 
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Figure 2: Phase portrait of the system ^^)-(^^) corresponding to the axisymmetric 
Bianchi type I cosmology with Am > 0. Note that W denotes a line of non-isolated 
equilibrium points. The dashed line represents the exact x = constant solution (l^j- 
l^^. See caption to Fig. 1 for notation. Arrows on the trajectories in this three- 
dimensional phase space have been suppressed since the direction of increasing time 
is clear using this notation. Grey lines represent typical trajectories found within 
the two-dimensional invariant sets, dashed black lines are those trajectories along the 
intersection of the invariant sets, and solid black lines are typical trajectories within 
the full three-dimensional phase space. 
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This implies that 



dz _9 / _ 16 
rfe "iH""" 27 



(36) 



and Eq. (|36D can be integrated exphcitly in terms of 0-time. 

In the following Section we determine the effects of a negative cosmological con- 
stant. This allows a direct comparison to be made with the Am > models considered 
above. 



IV Negative Cosmological Constant 

In the case where Am < 0, the generalized Friedmann constraint equation (|TUp implies 
that 

-Am> + QN^ > 0. (37) 

We may therefore normalize by employing the quantity y/ip'^ — Am- Defining the new 
variables 



u 



V^' - Am 

Vi^'' - Am 
6N^ 



w = 



V^2-Am' 

where < {u"^, v"^, w} < 1, and the new time variable 

d Id 
dE - Am d6 



(38) 
(39) 
(40) 



(41) 



implies that Eqs. (p!2|)-([l^) become 



dE 
dv 

dE 

dw 

dE 



u — w 



V3 + 



uv 



2v (^1 — u'^ — uij — VSu (l — V 



uv 

2 



(42) 
(43) 
(44) 



The phase space is bounded by the sets v = ±1 and w = 1 — u^, where the latter 
corresponds to a zero axion field. The dynamics is determined by the fact that the 
right-hand side of Eq. (|4^ ) is positive definite so that m is a monotonically increasing 
function. 
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IVa Isotropic Model for Am < 



In the invariant set w = 0, corresponding to the isotropic FRW model (/3 = 0), the 
system (H^)-(H3) reduces to the following two-dimensional system: 



i (l _ (2n2 - 1 - v^. 



(45) 
(46) 



dl 

dv 

dS 2 

The lines = 1 and f ^ = 1 are invariant sets, containing four equilibrium points 
Su,v These points are all saddles and are located at the intersections of the lines. 
Their eigenvalues are given by 

^1,1 : Ai = V3-l, A2 = -2(^3 + 1) (47) 

: Ai = 2(V3 - 1), As = -(V3 + 1) (48) 

^1,-1 : Ai = -2(V3-1), A2 = (V3 + 1) (49) 

Ai = -(V3-1), A2 = 2(V3 + 1). (50) 



5-1,-1 

The remaining two equilibrium points and their eigenvalues are 

1 \ - 1 - 10 



R : (u-,V- 



^'"7!j ^^"71' ^'"v^ 



A : {u+,v+) = 



A - ^ A - ^° 



(51) 
(52) 



Consequently, i? is a source and yl is a sink. Fig. 3 depicts the phase plane of the 
system 



IVb Anisotropic Model for Am < 

In the full system (^2]) -(^4]) with a non-trivial shear parameter, the equilibrium points 
and their respective eigenvalues are: 



W 



± 



V = ±1, u + w = 1; 
(Ai,A2,A3) = lo,V3 



1 



R : v = ■i=,u = —l,w = 0; 

v3 

(Ai,A2,A3) = -^(l,2,10) 
1 

A: V = —;=,u = l,w = 0: 

Vs 



(Ai, A2, A3) 



1,2,10). 



-2V3 



u± 



(53) 



(54) 



(55) 



V 




u 



Figure 3: Phase portrait of the system ^^-jj^) corresponding to the isotropic FRW 
model with Am < 0. See caption to Fig. 1. 

The saddle points S^i^-i in subsection 4.1 are the endpoints to the hne W~ . This 
hne represents early-time attracting solutions for — 1/\/3 < u < 1/v^ and saddles 
otherwise. The saddle points S±i^i are the endpoints to the hne . This corresponds 
to late-time attracting solutions for — 1/-\/3 < u < l/\/?> and saddles otherwise. 
Hence, there are two early-time attractors given by the point R and the line W~ 
for — 1/-\/3 < u < l/-\/3- There are also two late-time attractors corresponding to 
the point A and the line for — 1/-\/3 < u < l/-\/3. Fig. 4 depicts the three- 
dimensional phase space and table |^ lists each equilibrium set and its stability. 



Equilibrium Point 


Stability 


W-: Eq. (H) 


Repellor (Source) for ^ < m < 
Saddle otherwise 


W+: Eq. (H) 


Attractor (Sink) lox ^ < u < ^ 
Saddle otherwise 


R: Eq. (H) 


Repellor (Source) 


A. Eq. (H) 


Attractor (Sink) 



Table 2: Equilibrium sets for anisotropic model with Am < 0, and their stability (the 
equations where each set is defined is also listed). 

This concludes the derivation of the phase portraits for the spatially flat and 
homogeneous cosmologies derived from Eq. (^). We proceed in the following Section 
to discuss their properties. 
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w 
A 




Figure 4: Phase portrait of the system (^^-(^^ corresponding to the axisymmet- 
ric Bianchi type I model with Am < 0. Note that denote lines of non-isolated 
equilibrium points. See captions to Figs. 1 and 2. 
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V Interpretation of the Phase Portraits 



The dynamics of the isotropic cosmology described by the system (Pq)- (|2^) is of 
interest from a mathematical point of view due to the existence of the quasi-periodic 
behaviour. The orbits are future asymptotic to a heteroclinic cycle, consisting of 
the two saddle equilibrium points 5*1 and 5*2 and the single (boundary) orbits in the 
invariant sets z = and 1 — x'^ — z = joining Si and 5*2 (see Fig. 1). The former 
set corresponds to the zero Am solution (given by Eqs. ([15| ) with q = 0) and the 
latter to the solution with constant axion field (see Eq. (|3TD); to our knowledge this 
exact solution was not previously known. In a given 'cycle', an orbit spends a long 
time close to 5*1 and then moves quickly to 5*2 shadowing the orbit in the invariant 
set z = 0. It is then again quasi- stationary and remains close to the equilibrium 
point 5*2 before quickly moving back to 5*1 shadowing the orbit in the invariant set 
l — x"^ — z = 0. We stress that the motion is noi periodic, and on each successive cycle 
a given orbit spends more and more time in the neighbourhood of the equilibrium 
points 5*1 and S2. 

In Fig. 1, the exact solution corresponding to the equilibrium point F is a power- 
law solution: 

a = a,(-t)^/^ 



V 15Am, x2 

$ = 0, (56) 

where t is defined over the range — oo<t<Obya suitable choice of an integration 
constant. This new solution represents a cosmology that collapses monotonically to 
zero volume at t = 0. The curvature and coupling are both singular at this point. 
The universe is initially in a weak coupling regime, since Ges —>■ a.s t —00, and 
the effective energy density of the axion field also vanishes in this limit. 

All orbits in Fig. 1 begin at F. The cyclical nature of these orbits can be phys- 
ically understood by reinterpreting the axion field in terms of a membrane. The 
homogeneity of the axion field, a = <j(t), implies that the two-form potential, 5^,^, 
must be independent of cosmic time, and this in turn implies that its field strength 
must be proportional to the volume form of the three-space. If the topology of the 
spatial sections is given by a three-torus, x x S^, the behaviour of the axion 
field is dynamically equivalent to that of a membrane that has been wrapped around 
this torus The collapse is resisted by this membrane and the universe under- 

goes a bounce. As the volume increases, however, the influence of the membrane 
is diminished, because the energy density of the axion is rapidly redshifted away. 
Consequently, the cosmo logical constant becomes important. 
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The subsequent effect of the cosmological constant can be determined by view- 
ing Eq. in terms of a Brans-Dicke action, where the dilaton-graviton couphng 
parameter is given by = — 1 |Q. The Brans-Dicke FRW models containing only a 
cosmological constant in the matter sector have been discussed previously by Barrow 
and Maeda, but their solutions only apply for o; > —5/6 J^, The behaviour of 
the general solution for uj < —5/6 is different and can be established by performing a 
conformal transformation to a frame where the dilaton field is minimally coupled to 
gravity. In such a frame the term containing Am niay be viewed as an exponential, 
self-interaction potential for the dilaton, where the exponent is uniquely determined 
by the value of u ^^T^. When u > —5/6, the late-time attractor is a scaling solution, 
where the kinetic and potential energies of the dilaton field redshift in direct propor- 



tion For UJ < —5/6, however, the potential is so steep that the dilaton effectively 
becomes massless Thus, the late-time attractor when uj = —1 corresponds to 
the solution (|I6|) where hi = 1/3. 



Further insight may be gained by defining new variables in the reduced action (^): 

X = 4a - $ 

7 = $ — 6q; 

i= [ dte^. (57) 



In the case where = (t = 0, Eq. (§) reduces to 



S= dt 




(5^ 



The momentum conjugate to the variable x is constant, i.e., dx/dt = C, and the field 
equation for 7 is a Liouville equation: 

d^l 

^^Amc-. (59) 

The general solution to Eq. (|59D satisfying the Hamiltonian constraint can be found. 
When C > 0, it can be shown that 7 oc y/SCH in the late-time limit. Since the 
Hubble parameter is given by 



a = — C + ^ , (60) 




the late-time attractor corresponds to the collapsing solution in Eq. (p^ . 

In effect, therefore, the cosmological constant resists the expansion and ultimately 
causes the universe to recollapse and asymptotically approach the saddle point Si. On 
the other hand, the collapse causes the axion field to become relevant once more and 
a further bounce ensues. The process is then repeated with the universe undergoing a 
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series of bounces. The orbits move progressively closer towards the two saddles, 5*1,2, 
and spend increasingly more time near to these points. This behaviour is related to 
the fact that the kinetic energy of the shifted dilaton field increases monotonically 
with time, since Eq. (@) implies that (/3 > 0. 

When shear is included 7^ 0), F still represents the only source in the system. 
The orbits follow cyclical trajectories in the neighbourhood of the invariant set y = 
and they spiral outwards monotonically, since Eq. (^) imphes that dy/dQ > 0. After 
a finite (but arbitrarily large) number of cycles the kinetic energy associated with the 
shear parameter, j3, begins to dominate the axion and cosmological constant. The 
orbits then asymptote to the power-law solutions (p!6D. All orbits in the full three- 
dimensional phase space actually spiral outwards around the orbit represented by 
the dashed line in Fig. 2 which corresponds to the exact solution (|3^-(p6D with 
X = constant. In terms of cosmic time, t, this exact solution satisfies 



and 



16,2 48,^,2 

Am W H 

2V 13 



a = — 09, 







whence from Eqs. (|5|)- ([l0|) we obtain 

where k^p is an integration constant. Defining 



(61) 
(62) 

(63) 



g = 



(64) 



simplifies Eq. ( |63D to 

gg = —1 (65) 

and Eq. ( p^ can be integrated exactly to obtain A second integration then 
yields (f in terms of the Inverse Error function, so that in principle we can obtain the 
scale factor function of time, t, from Eq. (|62D . 

This cyclical behaviour does not arise if Am < (see Fig. 3). The equilibrium 
points A and R represent the power-law solutions: 



$* - In 




(J* 



(66) 
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where the + sign corresponds to the point R and the — sign to the point A. Initially 
the universe is collapsing and the axion field induces a bounce, but this field can 
not dominate the dynamics again once the volume of the universe has increased 
sufficiently. 

Fig. 4 depicts the axisymmetric Bianchi type I model when Am < 0. In this 
phase space, for —l/^/3 < u < 1/v^ the line W~ represents the negative branch of 
the solution (|T6]) for E (—1/3,1/3). Likewise, for —l/^/S < u < l/\/3 the line 
represents the "+" solution in Eq. (p!6D for G (—1/3,1/3). The four saddle 
points Su^v correspond to the power-law solutions (0) with = ±1/^/3. From Fig. 
4 we see that generically trajectories asymptote away from either the line W~ or 
the point R and move towards the expanding power-law solutions or A. Hence, 
the cosmological constant is important in determining both the early- and late-time 
dynamics. Since u is monotonically increasing (see (^21) ) we note that the occurrence 
of a bounce in these cosmological models is a typical feature. 



VI Discussion 

In this paper we have presented a qualitative analysis of spatially ffat FRW and 
Bianchi type I cosmologies containing non-trivial dilaton and axion fields with a 
cosmological constant in the matter sector of the theory. The action we considered 
reduces to the string effective action when the cosmological constant vanishes. A 
complete stability analysis was performed in all cases by ffnding variables that led to 
a compactiffcation of the phase space. We found that a cosmological constant has a 
significant effect on the dynamics of the string cosmologies (|T^. 

One of the more interesting mathematical features of the models we have con- 
sidered is the existence of quasi-periodic behaviour. This occurs in the isotropic 
cosmologies, where the orbits are future asymptotic to a heteroclinic cycle (see Fig. 
1). The solutions interpolate between the saddles 5*1 and 5*2 corresponding to the 
power-law models (|1^) with = l/y/S. It would be interesting to consider the 
implications of this behaviour for the pre-big bang infiationary scenario |T^. We 
note that the phase portrait depicted in Fig. 1 is similar to that of Fig. 1(e) in 



21| that describes the locally rotationally symmetric submanifold of the stationary 
Bianchi type I perfect fiuid models in general relativity, although in this latter case 
the independent variable is space-like. 

The general Bianchi cosmology, where the shear matrix is given by 

Pab = diag [p+ + V^P- , P+ - VSP. , -2/3+] (67) 

can be analysed directly by defining the variable N in Eqs. ( p!7| ) and (|20| ) via 
A^^ = J2i=±Nf. Orbits in the full phase space of Fig. 2 with non-trivial shear 
term (represented by the variable y) are repelled from the source F. The variable y 
increases monotonically and the orbits spiral around the exact solution given by Eqs. 
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(|61|)-(|63|), as represented by the dashed hne in Fig. 2. (See also Fig. 1(f) and the Ap- 
pendix in ||2l|). This imphes that solutions are asymptotic in the past to the solution 
given by Eq. ([561) . At early times the orbits 'shadow' the orbits in the invariant set 
y = and undertake cycles between the saddles (in three-dimensional phase space) 
on the equilibrium set W close to Si and 5*2. These saddles on W may be interpreted 



as Kasner-like solutions Note that y = at Si and 5*2, however, and there 

is no shear term in these cases. The orbits thus experience a finite number of cycles 
in which the solutions interpolate between different Kasner-like states. The orbits 
eventually asymptote towards a source on the line W. 

This is perhaps reminiscent of the mixmaster behaviour that occurs in the Bianchi 



type VIII and IX cosmologies [22, 23]. These are the most general models in the 



Bianchi class A of spatially homogeneous universes In these models, Taub orbits 
joining equilibrium points of the Kasner set K lead to the existence of infinite hetero- 
clinic sequences which approximate the past asymptotic behaviour of generic orbits. 
(These heteroclinic sequences are defined by a map of K onto itself). Mixmaster 
oscillations also occur in less general (i.e. lower-dimensional) Bianchi models with 
a magnetic field p5[ or Yang-Mills fields [^. It is interesting to note in the string 
context that mixmaster behaviour also occurs in scalar-tensor theories of gravity in 
general and in Brans-Dicke theory in particular p7 |. 

This analogy is only suggestive. We note that if a non-zero central charge deficit 
is included, the quasi-periodic behaviour in the full (higher-dimensional) phase space 
does indeed persist Unlike the mixmaster oscillations, however, the orbits in 

Fig. 2 eventually spiral away from y = 0, although there are orbits that experience a 
finite but arbitrarily large number of oscillations. However, it would be interesting to 
further explore any correspondence with possible mixmaster behaviour, particularly 
by including addtional anisotropic or matter degrees of freedom. 

Some of the dynamics discussed in this paper is also relevant to higher-dimensional 
cosmological models. Kaluza-Klein compactification of ten-dimensional supergravity 
theories onto an isotropic six-torus of radius e'^ introduces an additional modulus 
field into the effective four-dimensional action Integration over the spatial vari- 
ables for a spatially flat FRW model then leads to an action that is formally identical 
to that of Eq. (H) when we specify Am = 0. In this sense, therefore, the action (H) can 
be recast into a higher-dimensional context, where the shear term jS plays the role of 
the modulus fleld and Am is interpreted as a cosmological constant that is introduced 
after compactiflcation. 

More generally, type II supergravity theories contain Ramond-Ramond form- 
flelds that do not couple directly to the dilaton fleld in the string frame 0. Under 
dimensional reduction, these flelds give rise to terms in the effective action of the form 

exp(c/?), where Q and c are constants i.e., Ramond-Ramond charges give rise 
to exponential potentials for the modulus fleld rather than a simple constant term 
such as that considered in this work. However, from the analysis in section 3.1, there 
will be string solutions containing Ramond-Ramond flelds that asymptote towards 
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solutions with /? = constant {y = in Fig. 2), in which case Q^exp(c/5) is effectively 
constant. It might then be expected that the heterochnic cycle that occurs in the 
invariant set y = (see Fig. 1) will play an important role in describing the dynamics 
of these string cosmologies. 
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